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Abstract 

A behaviour of molecular cluster with Lennard- Jones potential of inter- 
actions as Hamiltonian system is studied by computer simulation (molec- 
ular dynamics method). It is shown that complex periodic oscillations 
of the cluster as a whole are possible. This is in accordance with KAM 
theorem. 

1 Introduction 

We study the behaviour of the dynamical system with Lennard - Jones interac- 
tion by computer simulation. The main physical motivation is concerned with 
long - time dynamical evolution in nonlinear Hamiltonian systems perturbed by 
random fluctuations caused by rounding-off errors. 

The goal of works [1-4] was the investigation of a dynamical foundations of 
classical statistical mechanics. It were considered : 

- a chain of iV nonlinearly coupled particles whose Hamiltonian is given by [3] 
( Fermi - Pasta - Ulam f3 - model ) 
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where <pi are the displacement with respect to stable positions, 
a chain of N nonlinearly coupled particles with Hamiltonian [3] 
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- the system of N particles arranged on a square lattice with nonlinear 
interaction through a Lennard - Jones potential [1] 
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e is a depth of the potential well; 
a = r 2- 1 / 6 , where V(r ) = 0, in [1-2] were taken Ai = 1, A 2 = 1 . 
Boundary conditions have been chosen periodic in normal mode coordinates [3] 

4>i = 4>i+N, 

or square lattice was surrounded by a border of fixed particles [1-2] . 

In all cited works it have been proposed that particles oscillate about their 
stable positions. The main results there obtained may be summarized as fol- 
lows: 

- at low energy there exist ordered motions, normal modes appear to be uncou- 
pled [2]; 

- corrcspondengly to a certain value of the average - time kinetic energy per 
particle relaxation toward the equilibrium is strongly slowed down [4]; 

- a time of relaxation r r (e) ( e is a full energy per particle in normal modes 
terms ) is compatible with Nekhoroshev - like low [4] with model dependent 
parameter 6 
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At high energy the transition between ordered - disordered ( weakly - strongly 
chaotic dynamics ) occurs but there are no exist thresholds between two dynam- 
ical regimes. In all cited papers the choice of an indicator of system's state was 
found to be insufficient. The use of relaxation time T r (e) is initial conditions 
dependent [4], spectral entropy r](t) is different from zero when t — > oo [3-4], Lay- 
pounov characreristic exponent A(e) is initial conditions independent but gives 
model - dependent transition and poor convergence when number of particles 
N is large. 

The main physical conclusion is the following: " there is a nonvanishing mea- 
sure of initial conditions in phase space that can originate long living metastable 
states even in the thermodynamical limit " [4]. 

It has been studied Hamiltonian dynamics subjected to rounding - off errors 
fluctuations. The tests concerned the Hamiltonian dynamics were the following: 
negligibly small value of Ae/e ~ 10~ 5 — 10~ 6 and time reversible dynamics for 
i--t [4]. 

Based on cited results it may be concluded that: 

- with the proviso that particles oscillate close to positions where their potential 
energy has a minimum, provided that number of the positions equal to number of 
particles , when full enrgy e > e c ( e c is model and initial conditions dependent) 
a number of oscillations with continuous random spectrum may increase on 
account of decrease in oscillations with discrete spectrum; 

- apparently in general way a characteristic time (e) for breaking of discrete 
spectrum is not existent; 

- a finite characteristic time rl 2 \e) within which one would to determine a 
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random nature of oscillations in system is not existent; 

- there is a fundamental problem of long - time following of dynamics in computer 
simulations because of rounding - off errors. This are nearly tangent to problem 
of modelling of dynamical evolution. 

The starting positions of the work based on assumption that weak changes of 
oscillations modes in dynamical system with large number of degrees of freedom 
may give rise to complex behavior of whole system when to abandon any restric- 
tion in motions and in particular one could to observe the dynamics without 
any changes in oscillations modes during long - time computer simulations with 
unbreaking correlations of the motions in a whole system despite of rounding - 
off errors. 



2 Model description 

Computer simulation runs have been performed for N particles Lennard -Jones 
system without any boundary or the like conditions. The equations of motion 
have been computed with undimensional time - step At — 0.05 — 0.13, where 
real time - step 
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At real = At ^ — J (5) 
The undimensional average - time kinetic energy per particle is defined as follow 
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Computer simulations were carried out without any cut off range of r. 
Initial conditions are of particular significance in present work. The starting 
points of particles are setting up as follows 

Xl = aAj sin (^^j , Vi = 4? cos (j^j (7) 

A = 6.752/4.53, N t = 22 when j = 4,i = 1...22, and N, = 32 when j = 5, 
i = 1...32, and N, = 44 when j = 6,i = 1...44, a = 1.013, A 1 = 4.870, A 2 = 1. 

At the first stage of formation of the cluster next scaling way was used. After 
time t = 200fcAt, k = 1...5, over the particles velocities were set v X i — 0, v y i = 
for i — 1...N. After some time t\ the scaling of velocities was abandoned and 
the cluster became as isolated but subjected to rounding - off errors. So the 
dynamical system is Hamiltonian or near Hamiltonian as fluctuations of full 
energy Ae/e were of the order of 10^ 4 — 10~ 5 and average - time (Ae) was equal 
to zero. 
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3 Results description 



The graphical displays of the cluster during program runs is depicted in Fig.l. 
The cluster oscillates as whole with period t* ~ 2000 and is unbreaking by 
rounding - off errors ( in this case the time of the unbreaking correlations of 
motion of cluster as a whole was t r = 5 10 4 , time step At = 0.13). The 
dependence of time of an unbreaking correlations vs its full energy per particle 
is reported in Fig. 2 ( ti = 1100 in this case ). In Fig. 3 is shown the value of 
average kinetic energy per particle. It is displayed time - average kinetic energy 
(ekin) though all particles of the cluster has a complex periodic oscillations about 
some state ( Fig. 4). 

The cluster doesn't lose this complex periodic oscillation when time - step of 
integration was 2 times increased and doesn't lose its symmetry when program 
runs ( t r = 10 5 ,time - step At = 0.13) in the case when full energy per particle 
is about 0.01 (in undimensional units). 

The number of particles with similar symmetrical behavior is no less then 
two - fold. 

4 Discussion 

In accordance with KAM theorem such cluster with complex periodic oscilla- 
tions is possible. 

What is to bring out the existence in computer simulation such cluster? 

- the cluster has sensitivity to initial conditions and to parameters Ai , A 2 of 
potential of interaction. All of this are characteristic for chaotic system. 

- the cluster brings out the unbreaking collective correlations. The average - 
time kinetic energy equal to 1-3 K ( parameters e, a were taken for argon ). This 
is the interval of an energy where convergence of normal modes in questionable 
in accordance with [1-2]. 

- collective correlations of the motion increase the stability of the cluster as a 
whole despite of rounding - off errors. 

Is the cluster artificial or has its physical reality ? This is out of the purpose 
of the present work. The dynamics of the cluster is near Hamiltonian with 
(e) ~ const. 

We would like to underline that symmetry in initial position of the particles 
brings complex periodic oscillations in this case. For the cluster with parameters 
A = 6.752/4.53, Ni = 20 when j = 4,i = 1... 20,7V* = 30 when j = 5,i = 
1...30, Ni = 40 when j = 6,i = 1...40, a = 1.000, A x = 4.870, A 2 = 1. one could 
to observe spontaneous turning of symmetry axis induced by rounding - off 
errors. Spontaneous breaking or turning of symmetry of any kind in Hamiltonian 
dynamics is impossible. 

Some words about description parameters of the cluster where one could 
observe contraction - extention of the cluster as a whole. When time - space scale 
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is large then one needs only 4 parameters to describe contraction - extention 
of the cluster as a whole. For molecular space and time scale one needs much 
more number of degrees of desription's parameters ( but less then 4iV). 

There is main question from the work. We would like to formulate this ques- 
tion as follow: What number of degrees of freedom, full energy and potential of 
interaction did provide time of an unbrcaking correlations in a cluster subjected 
to rounding - off errors for an infinity or this time is always finite ? 
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